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Abstract. Rough estimates and upper limits on current 
and particle densities form the basis of most of the canon- 
ical pulsar models. Whereas the surface of the rotating 
neutron star is capable of supplying sufficient charges to 
provide a current that, given the polar cap potential, 
could easily fuel the observed energy loss processes, ob- 
servational and theoretical constraints provide strict up- 
per limits to the charge densities. The space charge of a 
current consisting solely of particles having only one sign 
creates a compensating potential that will make the max- 
imum current dependent on potential and distance. In the 
non-relativistic case this fact is expressed in the famil- 
iar Child-Langmuir law. Its relativistic generalization and 
subsequent application to the inner pulsar magnetosphere 
provides clear limits on the strength and radial extension 
of charged currents originating on the polar cap. Violent 
Pierce-type oscillations set in, if one attempts to inject 
more current than the space charge limit into a given vol- 
ume. These considerations apply wherever there is a sig- 
nificant amount of charged current flow, in particular in 
the gap regions. There they can be used to derive limits 
on the size of such gaps and their stability. 



1. Introduction 

All pulsar models require charges to be drawn from the 
pulsar surface, accelerated, and possibly multiplied by pair 
production, as being the medium that will be responsible 
for the great variety of observational features currently 
known from many pulsars. Common to all models is a 
rotating heavily magnetized (10^ — lO^^G) neutron star. 
As early as 1955 it was recognized by A. Deutseh that a 
sufficiently strongly magnetized fast rotator in a vacuum 
would be surrounded by an envelope (magnetosphere) of 
charged particles. Because of their fundamental impor- 
tance, the electrodynamical aspects of pulsar models have 
been of interest from the early days of pulsar research until 
now. 



In all cases, the number density of particles available 
is the decisive factor for the self-consistency and hence 
viability of any pulsar model. 

Starting from the surface of the neutron star, it is rec- 
ognized that its very composition plays a crucial role in 
determining the number charges emitted (Usov & Melrose 
1996). The controlling factor is the surface work func- 
tion, representing the potential barrier to be overcome 
by the charges that are to be emitted from the pulsar 
surface. Usually, the work function is of the order of the 
Fermi-energy of the solid lattice, which in itself depends 
on the electron density. Using a Thomas-Fermi-Dirac- 
Weizsacker approximation Abrahams & Shapiro (1992) 
calculated the cold (T = 0) surface density ppe to be 
about 2.9 • lO^kgm^'^. That result enabled us (Jessner et. 
al. 2001) to estimate the Fermi energy at the surface as 

2 ■ f pFe-{28-2) Y 
^''^P'^^ = ^^B^ ■ [ 56. ) 

which amounts to Ep = 4.17 • 10^ eV. For such a low work 
function, even field emission can supply enough charges to 
provide a relativistic flow of electrons at the co-rotation 
value given by the Goldreich- Julian-density (Jessner et. al. 
2001). But the estimated value for Ep should be taken as 
an upper limit, since the thermal X-ray emission indicates 
that pulsars have surface temperatures well above 3 • lO^K 
and are surrounded by a thin gaseous surface layer which 
would not entail any significant barrier for the extraction 
of electrons (Pavlov and Zavlin 1998). In that case the sur- 
face properties themselves will not influence the primary 
charge densities (free charge emission scenario). 

Recently (Jessner et. al 2001), we proceeded to calcu- 
late the Lorentz factors and pair production rates of elec- 
trons accelerated along field lines of an aligned dipole tak- 
ing due account of radiation losses due to inverse Comp- 
ton scattering on thermal photons (as was proposed by 
Sturner 1995) and curvature radiation. This was done un- 
der the simplifying assumption of a relativistic flow which 
resulted in constant shielding of the accelerating fleld by 
the primary charges below the pair production region. In 
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Fig. 1. Simple one dimensional flow of electrons between 
two electrodes. 

addition, the field and density enhancements described by 
Muslimov and Tsygan (1992) were also not implemented. 
To be complete we will now give an estimate of the initial 
(primary) current that can be drawn from a pulsar polar 
cap and on its stability by solving Poisson's equation with 
constant current density, but with a velocity dependent 
charge density. 

2. The Child-Langmuir Equation 

Space charge limitations play an important role in the con- 
duction of currents across otherwise empty space, their 
detailed study in first half of the 20**^ century provided 
the foundation of much of our electronic technology. In 
order to illustrate the principle, we will briefly review a 
simple classical 1-dimensional case: Take two large elec- 
trodes (Fig. 1) of equal area Ao separated by a distance x, 
one at zero potential and the other at U = Uq. We assume 
that we can supply an arbitrary amount of charge for the 
flow in the direction of x, so that the current will limit it- 
self as soon as the potential at the lower surface vanishes. 
The velocity of an electron, accelerated by a potential U 
is given by 



PiU) 



(2) 



In the non-relativistic case {U « 511 kV) this sim- 
plifies to 



v{U) 



2-qeU 



(3) 



If the potential U > 0, then a current / will flow from 
the lower to the upper electrode. This current is conserved 
along the upward (x) direction. Such a current moving in 
a vacuum obeys Poisson's law: 

I 

U = (4) 



In the non-relativistic case we replace cl3 with . / ^'^g^ 
integrate twice and solve for the current to obtain the 
canonical expressions for the maximum current Iq and the 
internal resistance Rq = U/Ia of a planar vacuum diode 
(Child, 1911, Langmuir 1913) that has a potential U w.r.t. 
the cathode. 



Io{U,x) = ^.[^eoA,U^/^x-' 
6 V nie 



The resistance follows as 
Ro{U,x) 



4^/2eo^01 



= [/" 



(5) 



(6) 



Hence a very important, but often overlooked conse- 
quence of Poisson's equation and the equation of motion 
is the fact, that a charged current requires a potential to 
flow through a vacuum. In the non-relativistic regime, the 
required potential increases as the 4/3 power of the longi- 
tudinal distance and the 2/3 power of the current. 

The solution to the fully relativistic equation 



U 



I ■ {qeU + TUeC^) 



is slightly more complicated, 



J((7, x) = ( 2^qlU^ + 2geC/mec2 

qex^ y 

-TTipC 2 sm 



(7) 



(8) 



dx"^ 



eocAoPiU) 



qeU + rrieC^ 



but the general dependence of the space charge lim- 
ited current on the inverse squared distance remains un- 
changed. In fig. 2 we give a quantitative example of the 
current that will flow between two electrodes of Icm^ 
crossection, spaced apart by one meter. In the relativis- 
tic regime, the current rises linearly with the applied 
potential, hence the potential required to drive a par- 
ticular current density over a distance x now increases 
as .T^. Although this example appears to be somewhat 
remote from the expected conditions at a pulsar's po- 
lar cap, one easily recognizes, that for a typical Pulsar 
(B = lO^^G, P = 0.5s), the relativistic Goldreich-Juhan 
current density approaches ^ qen{rns)c = 498 Acm^^. 
A potential of 10^ V is required to drive such a cur- 
rent over the distance of 1 meter. For a height of i.e. 
h — 100m we have exhausted the available potential dif- 
ference (~ lO^'^ V) on the polar cap! The square law de- 
pendence of the required potential makes it impossible to 
have a far reaching charged current of appreciable density 
within the pulsar magnetosphere. 

^ Here we used the general relativistic surface density n{r„s) 
as defined in section 3 
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Child - Langmuir Law 




Potential U (V) 



Fig. 2. Current flowing between two parallel Icm^ elec- 
trodes ( spacing 1 m) as a function of the potential dif- 
ference U (solid line). The dashed line shows the usual 

2 

non-rclativistic relationship which is invalid for U > . 




Fig. 3. Aligned magnetic rotator 



3. Space charge limits for an aligned relativistic 
pulsar magnetosphere 

For reasons of simplicity we will now use the model of 
an aligned rotator as described by Goldreich and Julian 
(1969). Fig. 3 outlines the geometry of the rotating neu- 
tron star field. Near the surface of the neutron star, the 
space-time metric is affected by the presence of the heavy 
compact rotator leading to the general relativistic effects 
of the Lense-Thirring dragging of field lines and a radial 
dependence of the volume elements. Muslimov and Tsy- 
gan (1992) introduced an auxiliary function to account for 



the general relativistic distortion of the field lines: 

3 r 



f{r) = -3 



In 1 



(9) 



2r 



here we used the compactness cmt 



2rM„ 



where 



2rM„3 
— 3 — 



Tg is the canonical Schwarzschild radius of the 



neutron star. A magnetic field line starting with a colat- 
cral angle (j)o at the pulsar surface is given by 



sin(0) 



/(I) 



(10) 



As an example we use the well known pulsar PSR 
B2021-I-51 from the catalogue of Taylor, Manchester and 
Lyne (1995) with a period of P = 0.529s and a surface 
magnetic field of Bq = 10^^'^^G which for simplicity we 
treat as if it were an aligned rotator. Then the colat- 
eral angle of the rim of the polar cap turns out to be 

J, ,,:„-l 



^W2) 



0.94 and the area of 
c„t)) = 8.5-104m2. 



Pc/(1 

the polar cap is Apoie = 27rr^^(l— cos 
The opening angle (j> and thus the surface area enclosed by 
the limiting field lines increases with increasing distance 
from the surface. Using the equation of field fines (10), we 
can subsequenty replace cos(^(r)) with 



C(r) = 



\ 



1 - 



/(I) 



/ 



t) (11) 



The Goldreich-Julian charge density compensates the in- 
duced electric field parallel to the magnetic field line and 
is classically given for the pole by no{r) ~ ^<>^^^a (^)'^ 
yielding no(r„s) = 8.5 • lO^^m"^ on the surface of the cen- 
tre of the polar cap. Following Beskin (1999) we can use 



fo{r) 



i-(^r 



(12) 



1 



to describe the G.R. correction of the co-rotation density 
n as a function of distance. Because of the higher electric 
field at the surface, the density has to increase by a factor 
of fo{rns) = 1-21 to a value of no(r„s) • /o(r-„s) = 1.03 • 
10^''m"'^ . fo{r) decays quickly with r, i.e. to 1.021 at 
r = 10r„s. Assuming relativistic velocities, the total polar 
cap current can roughly account for the observed spin- 
down energy loss of the pulsar: 



leap = n{rns)qecAo- 



1 



'■'crit Jo 



(3cos'^{e)-i)de (13) 



1.01 • lO^^A 
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Together with the polar cap potential Ucap = 
{cos^{(j)c„t) - 1) = 2.07-1013V we obtain an elec- 
trical power of IcapUcap = 2.09 • lO^^W which, for such a 
coarse model, agrees surprisingly well with the observed 
energy loss of 1032-9££g = 7.9 . 1025W. 

For a pulsar with an electron current flowing out of the 
polar cap, the radial part of Poisson's equation can now 
be written in spherical polar coordinates: 

1 <^ 2 rr M A\ 

r — U = (14) 

r-^dr dr eo27rr2,c/3(u)(l - C(r)) ^ ' 

+ ^no(r„.)/o(r)(^)'(3C(rf-l) 
eo \ r J 

The first term on the right hand side of the equation con- 
tains the charge density due to the current of electrons and 
the second term the equivalent co-rotation charge density. 
Here we have assumed a uniform current density across 
the cap and we neglect the transverse dependence of the 
potential (all other derivatives in Poission's equation are 
assumed to vanish). But the divergence of field lines as a 
function of distance is accounted for by C,{r). Furthermore, 
it is quite reasonable to assume that the electrons will 
leave the surface with thermal velocities. For T = 3 • lO^K 
their kinetic energy corresponds to an initial potential of 
Uo = Y' = 25.9V, and a velocity of 0.01c. Come what 
may, the boundary condition at the surface < U < Uq 
dictates that the flow starts in the non-relativistic regime. 
As current densities are conserved along the field lines, 
and charges are accelerated, the first (current) term in 
(14) will decrease on acceleration, leading to further in- 
crease of the accelerating potential (and vice versa). If we 
therefore have a space charge limited fiow with full shield- 
ing at the surface (n = n(r„s)) and /3 = 0.01, the charge 
density will decrease by a factor of ~ /? along the field lines 
if relativistic potentials arc encountered, a point already 
made by Fawley, Arons and Scharlcman (1977). 

4. Numerical Solutions of the Pulsar Poisson 
Equation 

A Buhrsch-Stoer algorithm (provided by MATHCAD) 
was employed with typically 2048 steps to solve (14) 
for a range of up to 6m above the pulsar surface where 
the boundary conditions were set to U (r„s ) = Uq and 
^ = OVm~^. The results were checked by increasing the 
resolution and sometimes with another algorithm (adapt- 
able Runge-Kutta method). 

4-.1. Strict space charge limitation at the surface 

With a constant current over the polar cap that is in 
equilibrium at the surface, fulfilling the above boundary 
conditions, we get a massive increase in particle energies 
with height (fig. 4). But in the same vein, the charge den- 
sity falls to about 1% of the surface value after about 
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Fig. 4. Lorentz factor reached by particles in a current 
that is limited by the space charge density being equal to 
the co-rotation density at the pulsar surface. 

^ = cy^ ^°2™° = 15mm in agreement with Fawley, Arons 

and Scharleman (1977). In this and all following cases we 
neglected energy losses of particles by inverse Compton 
scattering and curvature radiation. These losses will start 
to limit the attainable particle energies in the surface layer 
as soon as Lorentz factors of about a few 1000 are attained 
(Jessner et.al. 2001). 

4-2. Enforced current at about 0.8 times the relativistic 
co-rotational current 

The extremely low relativistic charge densities of 
0.01n(r„s) obtained from the steady-state solution make 
pulsar models unworkable, if they require a primary cur- 
rent of co-rotational density (Ruderman & Sutherland 
1975) to supply the energy for the magnetospheric pro- 
cesses. Following the reasoning in section 3 we now find 
that picapVcap is more than two orders of magnitude 
smaller than the observed energy loss. The high, un- 
shielded may aid pair creation, but unless the sec- 
ondary charges themselves are accelerated by an un- 
shielded field, they will not contribute to the energy bal- 
ance. 

It is therefore of interest to see what happens when 
the strict boundary condition i7(r„s) = Uq is relaxed by 
forcing a charge overdensity to be drawn from the surface 
yielding a relativistic current density of the order of the 
co-rotation density. Fig. 5 shows the strongly fiuctuating 
potential distribution that will be the result. The current 
breaks into slabs of ~ 10cm thickness until the particles 
have obtained enough energy so that their charge density 
becomes low enough for the accelerating field to take over, 
in this case after ~ 5m. Such oscillations are common in 
fiows that exceed the space charge limit locally and the 
stationary solutions of similar cases have been described 
by Mestel et. al. (Mestel et. al. 1985) and Shibata (1997). 
Haeff noted the time variability of the instability in 1939, 
Pierce described it in 1949 and Eilek presented (this vol- 
ume) a time dependent solution showing how the oscil- 
latory patterns evolve and propagate with time. Schop- 
per presented a 3-D PIC simulation (this volume) show- 
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Fig. 5. Effective potential along a central field line for an 
enforced current of about 0.8 of the relativistic co-rotation 
current Iq. 

ing how an initially steady particle beam breaks into the 
noted charged slabs. 

4-3. Enforced currents with an additional positively 
charged relativistic inflow 

The local charge ovcrdensity on the surface, required by 
strong relativistic currents, can be neutralized by an in- 
flowing current of the same magnitude, but with opposite 
sign of the charges. The simplest case would be a highly 
relativistic current Ir with ^m^c^ >> qeU . Then we can 
set /? = 1 for the return current and solve the following 
differential equation: 



dr dr 



-1 



I 



-Ir 



(15) 



-noirns)fo{r) 



ns 

r 



(3C(r) 



1) 



The positively charged current will not neutralize the 
space charge in an extended region above the polar cap, 
leading to strong acceleration, similar to that shown in fig. 

4. But, not unexpectedly, as soon as one relaxes the neu- 
trality condition at the surface, density oscillations will 
appear. For closely matched inflowing and outflowing cur- 
rents a second boundary with U = Q will exist higher 
up in the magnetosphere. Fig. 6 shows an example where 
the effects of a positively charged current of 0.9/o and a 
negative thermal current of I.S/q are calculated. Further 
simulations with different flowrates have shown that the 
location of the second boundary is very sensitive to the 
matching of the currents. A decrease of the outflow by 
one percent moves the second boundary outward by sev- 
eral hundred meters. Figs. 7 and 8 show the oscillations 
near the surface. Mildly relativistic potential fluctuations 
of a few lO^V are encountered with wavelengths of about 
40cm and the densities (Fig. 8) show positive values with 
narrow gaps. 

5. Conclusions 

Charged currents cannot propagate in a vacuum without 
an outer potential, for a fixed current density the required 
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Fig. 6. Lorentz factor of outflowing particles with 1.8 of 
the relativistic co-rotation current /q and a positive inflow 
of 0.9/0. Particles come to rest sX, h ^ 2200m! 




height (m) 



Fig. 7. Potential oscillations caused by outflowing parti- 
cles with 1.8 of the relativistic co-rotation current /q and 
a positive inflow of 0.9/o- 
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Fig. 8. Charge densities for outflowing particles with 1.8 
of the relativistic co-rotation current Iq and a positive 
inflow of 0.9/0. 



potential increases as the square of the distance. In the 
case of pulsars, the available potential on the polar cap 
allows for a relativistic electron current of co-rotation den- 
sity only up to a height of a few hundred metres. 
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Particles start from the neutron star with thermal ve- 
locities. An electron current, if limited by non-relativistic 
space charges at the surface, will encounter an outwardly 
increasing potential along the field lines. This accelerates 
the charges and decreases the available current density. 
Hence a thermal current of non-relativistic surface den- 
sity cannot compensate the parallel electric field once it 
leaves the surface. In a steady-state regime, relativistic 
currents from the polar cap can only be a small fraction 
0(0.01) of a co-rotational current. Standard pulsar mod- 
els assume a relativistic current of co-rotational density 
and become unworkable when the primary currents are as 
low as shown. The available power in the magnetospheric 
current system is insufficient if the current is space charge 
limited. Due to the unshielded Lorentz factors would 
become too high for radio emission at the observed heights 
and if pair production were to take place, another boost 
of at least a hundred in the multiplicity factor would be 
required to reach the frequently assumed 10'^~^no(r). We 
take this as evidence against a steady state description of 
the pulsar mechanism and a large scale charged current 
flow out to i.e. the light cylinder. 

Forcing a current that locally exceeds the co-rotation 
charge density leads to strong Pierce type oscillations. The 
presence of a relativistic inflow of opposite charges does 
not alter the picture significantly, mainly the outflowing 
density can be increased by the inflowing amount. A sec- 
ond boundary with U = (inner gap) can be formed for 
particular combinations of inflow and outflow, where the 
outflow exceeds the inflow. In these cases the current and 
the potentials fluctuate with wavelengths of a few 10cm 
in a region of a few metres above the surface. The density 
is strongly modulated, showing narrow gaps in the flow. 
Such fluctuations are inevitable whenever the density ex- 
ceeds the space charge limits. It may happen at a bound- 
ary as in our case, or whenever a relativistic current is 
strongly perturbed so that it may become non-relativistic 
in places and by that way exceed the space charge limit. 
Work by Eilek (this volume), Schopper (also this volume) 
, Shibata (this volume and 1995), Mestel et. al. (1985) 
docs indicate that the Pierce instability is of importance 
in the pulsar magnetosphere. 
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